Light reflection upon a movable mirror as a paradigm for continuous variable 

teleportation network 
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We present an optomechanical system as a paradigm of three-mode teleportation network. Quan- 
tum state transfer among optical and vibrational modes becomes possible by exploiting correlations 
established by radiation pressure. 
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I. INTRODUCTION 

Quantum teleportation is one of the most fascinating 
possibilities offered by quantum information processing 
jlj . In the standard protocol, Alice is to transfer an un- 
known quantum state to Bob using, as the sole resources, 
some previous shared entanglement (quantum channel) 
and a classical channel capable of communicating mea- 
surement results. Often the quantum channel for con- 
tinuous variable (CV) teleportation is realized by using 
two-mode squeezed (TMS) states which may mimic non- 
classical Einstein-Podolsky-Rosen (EPR) correlations . 
The concept of quantum channel can be generalized to 
more than two parties to form a teleportation network. 
In such a case it can be used for telecloning 3], or, again, 
for teleportation by distilling bipartite entanglement [^. 

As a consequence of a famous science fiction saga, a 
wish to involve macroscopic systems into teleportation 
protocols now hovers in people mind. Recently, we have 
gone along this way showing the possibility to teleport 
a quantum state onto the vibrational mode of a mov- 
able mirror p|. In doing so we have exploited an intense 
radiation field impinging on the movable mirror and re- 
flected into an intense elastic carrier together with anclas- 
tic sideband modes. The latter together with the mirror 
vibrational mode constitute a tripartite quantum system. 
Here we shall examine it as a teleportation network where 
we exploit the bipartite entanglement, and more strongly 
the EPR correlations, which can be distilled for certain 
groupings of two modes, simply tracing out the remaining 
mode or measuring it by heterodyne detection and com- 
municating the result. In both situations the distilled 
channel comes out as a (nonzero mean) Gaussian bipar- 
tite state. The model is presented in Section |n] Then, 
we report the general teleportation protocol through a 
Gaussian channel in Section ITTT1 Hence, we examine the 
two situations above in Sections IIVI and Ivl Section IVT1 is 
for conclusions. 



II. EFFECTIVE HAMILTONIAN AND SYSTEM 
DYNAMICS 

We consider a perfectly reflecting mirror and an intense 
quasi-monochromatic laser beam impinging on its surface 
(see Fig. ^| . The laser beam is linearly polarized along 
the mirror surface and focused in such a way as to excite 



Gaussian acoustic modes of the mirror. These modes de- 
scribe small elastic deformations of the mirror along the 
direction orthogonal to its surface and are characterized 
by a small waist, a large quality factor and a small effec- 
tive mass @. It is possible to adopt a single vibrational 
mode description limiting detection bandwidth to include 
a single mechanical resonance of frequency CI. In this de- 
scription the incident laser beam, with frequency ojq, is 
reflected into an elastic carrier mode, with the same fre- 
quency loq, and two additional weak anelastic sideband 
modes with frequencies loq ± Cl. The physical process 
is very similar to a stimulated Brillouin scattering, even 
though in this case the Stokes and anti-Stokes component 
are back-scattered by the acoustic wave at reflection, and 
the optomechanical coupling is provided by the radiation 
pressure. Treating classically the intense incident beam 
(and the carrier mode), the quantum system is composed 
by three interacting quantized bosonic modes, i.e. the vi- 
brational mode and the two sideband modes. In our de- 
scription, vibrational, Stokes and anti-Stokes modes are 
denoted as 0, 1, and 2 respectively. In general the k th 
mode (k — 0,1,2) is characterized by ladder operators 
cik, a\ and by quadrature operators X k = (afc + q\:)/V2, 
P k = (a k - a[)/iV2 {[X k ,P k ] = i). In HQ] we have de- 
rived an effective interaction Hamiltonian for that system 



H e ff = -ihx{a\a 



a[al) 



ih9(a 2 al - a\ao) , (1) 



where \ and 9 are couplings constants whose ratio r = 

1/2 

6 fx = [(^o + CI)/(luq — Cl)] ' > 1 only depends on the 
involved frequencies. The system dynamics is satisfacto- 
rily reproduced by the Hamiltonian of Eq. £Q) as long as 
the dissipative coupling of the mirror vibrational mode 
with its environment is negligible. This happens if the in- 
teraction time is much smaller than the relaxation time of 
the vibrational mode (which can be ~ 1 s ||) and there- 
fore means having a high-Q vibrational mode (typically 
CI ~ MHz). 

The dynamics can be easily studied in terms of 
the symmetrically ordered characteristic function $(/x), 

[i = (^ ,Mi>M2) 9], where [i k = (jfj*® + ijj>^ is 
a complex variable corresponding to the operator a k 
(k = 0,1,2). The relation between the density op- 
erator and the corresponding characteristic function is 



given by p i2 
D k (n k ) 



[ULo / <PvkDl(»k)MlJ), where 



exp{/z fe a| - nl&k} = exp{i\/2(^ Xk 



2 



Mfe Pk)}- Notice that the Fourier transform dehned 
as above through the displacement operators -Dfc(/ifc), 

creates the correspondence fXjP «-» Xk, — Mfc *~* Pk- 
From Eq. we can deduce the dynamical equation for 
t). If we assume the initial condition Q(p, t — 0) — 



exp 



-n\IM>\' ~ ELo iMfcl 2 / 2 
uum state for modes 1,2 and to a thermal state with 
mean thermal number of excitations n for mode 0, the 
total state at time t is Gaussian, with characteristic func- 
tion 

$(M,t) = exp[-Q |W)| 2 - Qil^il 2 - Q2IM2I 2 + Jo(/ii^ 2 
+ M1M2) + 2iO«oM2 + MSM2) + T 2 {fi 1 [i + M1M0)]) ( 2 ) 

where: Q = B + 1/2, Q x = .4 + 1/2, Q 2 = £ + 1/2, 
T = T,T X = V,T 2 =C and the coefficients .A, B, C, 
V, £ and JF, explicitly given in 0, depend upon n, r 
and the scaled time t' - 



corresponding to a vac- 



2 . Since the state of 



Eq. is a zero-mean Gaussian state, it can be fully 
described by its correlation matrix (CM) V, dehned by 
Virn = (A£ ; A£ m + Af m A&)/2 (l,m = 1,...,6), where 
A£i = £i — and £ denotes the vector of quadratures: 
£ = (Xq, Pq,Xi, Pi, X 2 , P 2 )- In fact, introducing the real 
vectors fa (k = 0, 1, 2), defined by 



P-k € 



(3) 



so that C 3 3 fj, = (hq, Hi,fi2) < — ► (fa, Pi, fa) = P G K 6 , 
and expressing <i>(/i, i) in terms of ft, we get from Eq. (J2J 

$02,i)=e-^ T , (4) 
where the CM V appears and it is explicitly given by 
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(5) 



Since we deal with a closed system of three interacting os- 
cillators its dynamics is periodic in t' with period 2it. The 
separability properties of the tripartite system p 012 and 
those of the bipartite reduced systems = trk{po\ 2 ) 
(k = 0,1,2), have been already studied in Q. Here we 
briefly recall these properties referring to ^3 f° r the defi- 
nition of entanglement classes of a tripartite CV Gaussian 
state. The state poi 2 is almost everywhere fully entan- 
gled (class 1) except for isolated times t' — 2mir (m G N) 
when it is fully separable (class 5) and t' = (2m + 1)tt 
when it is one- mode biseparable (class 2). In the latter 
case, it is equal to the tensor product of a TMS state for 
the optical modes and of a thermal state for the vibra- 
tional mode. 

The inseparability between one of the three modes and 
the other two parties can be exploited to implement a 



telecloning protocol 0. When we trace out one of the 
three modes, we can distill bipartite entanglement only 
between modes 1 and 2 or between modes 1 and 0. In 
the first case (1 and 2), the entanglement between the 
optical modes exists at all times (except t' = 2mn), 
and it is extremely robust with respect to temperature. 
Such modes show also robust EPR correlations which 
are temperature-independent at t' = (2m + 1)tt when the 
state is a TMS state 0- In the second case, the entan- 
glement between the optical Stokes mode and the mirror 
vibrational mode exists in two limited time intervals, just 
after t' = and just before t' = 2tt, and they become 
narrower by increasing the temperature. In [5j the sec- 
ond time interval has been exploited to realize quantum 
teleportation from an optical to the mirror vibrational 
mode. 



III. CV TELEPORTATION THROUGH AN EPR 
CHANNEL 

Suppose to have a quantum teleportation network, i.e. 
a quantum channel given by a truly multipartite entan- 
gled state shared among N parties. Such a situation has 
been studied in |4j where, from a particular N— partite 
entangled state, a bipartite entanglement between any 
two of the N parties can be distilled to enable quantum 
teleportation. Clearly the distillability of bipartite entan- 
glement from the total channel is a necessary condition 
to make the teleportation network really quantum (i.e. 
not reproducible by any local classical means). 

Here we consider a distilled channel consisting of a bi- 
partite Gaussian state with a known drift. That happens 
for example when, from a iV-partite zero-mean Gaussian 
state, we trace out N — 2 modes, or when we measure 
them by heterodyne detection (communicating the re- 
sults to Bob through a classical channel). Suppose that 
such a distilled channel is shared in a network by Al- 
ice (mode i) and Bob (mode j) (for simplicity consider 
N = 3 and see Fig.0). If their channel is separable then 
they cannot perform quantum teleportation but, if it is 
entangled, then they can try to perform it by exploiting 
possible EPR correlations of their channel. We define 
"EPR+" or "EPR—" correlations in the following way 

EPR± (AXl) + (AP^) < 2, (6) 

where X ± = (X l ± Xj) and P ± = (P, ± Pj). De- 
pending on the supposed EPR correlations, Alice and 
Bob can implement an appropriate teleportation proto- 
col. We treat both cases in a compact form adopting 
the phase-space approach of The distilled channel 
is described by a Gaussian Wigner function W ch (at , ctj ) 
with a known drift and supposing to possess EPR± cor- 
relations according to definition JBJ. The unknown in- 
put state at Alice's station is described by W m ("f) so 
that the total state before the beam-splitter is given 
by W(j,ai,aj) = W in (j)W ch (a i ,a J ). After the beam 



3 



splitter (Fig. 0) we have 



v + — V- 

~7T 



w 



ch 



V2 



(7) 

where v± = (aj ±7)/ V% are the complex amplitudes 
of the output modes ± having quadratures x± = (Xi ± 
X in )l\[2 and p± = (Pi ± Pi n )/\/2. Depending on the 
(supposed) EPR correlations of the channel, Alice mea- 
sures an appropriate pair of quadratures of the output 
modes, namely EPR± implies measuring (x±,p^_). In 
terms of complex amplitudes, Alice's measurement re- 
sults (x±,p^) are expressed by the variable 7j_ with prob- 
ability density P("f'±) where 



1 

7+ 



V2(v 



(R) 



and j'_ = y2(i 



,(R)x 



Alice's detection collapses the total state {7J) whence 
Bob's conditioned state results 



W{ aj \i ± ) = 



1 



#lW*{>i)W A [±{il->f),a i 



(8) 

State (jSJ) is exactly the state at Bob's station when Bob 
receives from Alice her result -f' ± through a classical chan- 
nel. At this point Bob completes the teleportation pro- 
tocol by a suitable displacement of his mode, namely 



7± - 



(9) 



where the first quantity 7^ balances Alice's detection 
while the second one S± (to be specified) balances the 
drift of the channel and depends on the type of EPR cor- 
relations, too. The final state at Bob's station is then 
given by W(ctj — j± — 6±\j'±). In order to get the tele- 
portation fidelity, we compute the mean state teleported 
to Bob by averaging over all possible results 7^ 

W° ± u \a 3] 8 ± ) = J dY ± P(Y ± )W(«j - 7 ' ± - 5 ± |V ± ) 

= J d 2 1 K ± (a J , r ,S±)W m ( 1 ), (10) 

where 

K±{a j)T ,S ± ) = J d 2 7 ' ± W ch [±( 1 £-'y*),a j -Y±-8±] 

(11) 

depends also on the parameter 5±. The input-output 
relation (|10f> can be greatly simplified if we introduce 
the symmetrically ordered characteristic functions given 
by the Fourier transform of the Wigner functions, i.e. 
$(A) = J d 2 ae Xa "- x ' a W(a) where A = \W + ikW is 
the conjugate variable of the complex amplitude a. In 
fact, applying the Fourier transform to QlUp. we get: 



_ Ail -A*<5± ffydi 



$ cft ( T A*,A)$ m (A). (12) 



where ^ m (/i in ), <& out (iij) and 3> ch (fi t , [ij) are evaluated 
according to the substitutions shown in Eq. I|12|) . If we 



consider a pure input state then the fidelity takes the 
form 



d z \ $ m (A) [$ c/l (TA*,A)]*e 



-\5* ± +\"8± 



F±(6±) = - 

7T 

(13) 

where we have used Eq. i|12|l . In particular we con- 
sider a Gaussian pure input state and it can always 
be chosen with zero-mean, because of the invariance 
of F±(S±) with respect to the amplitude of the input 
state. Using the real variables of Eq. (J2J), the char- 
acteristic functions of the input and the channel are 
given respectively by $ m (# in ) = exp(-/Z in F m ^) and 
= exp[-(ft,/Z J 0V r '*(ft,/2j) T + D ch ( f I t , f l j ) T }, 
where the CMs V m and V ch are real positive matrices, 
while the drift vector D ch has pure imaginary elements. 
In particular we can write: 



A C 
C T B 



D cn = 2i (dt ,d 2 ,d 3 , d 4 ) 

with dk E M. Inserting $ in (/Zi„) and <& ch (fli,fi ) in 
Eq. lO, we get: 



F ± («5 ± ) = -=_=exp(-Q ± ) (14) 
v /det(E ± ) 



where 



E± = 2V m + RAR + B±RC± C T R T (15) 
/ 1 \ 



R = 



-1 



and 



}±^D ± (E ± )- 1 (D ± ) T >0 



D-t 



d 3 ,-S^±d 2 



di). 



(16) 



(17) 
(18) 



In Eq. (|14() the positive matrix E± depends only upon the 
CMs V m , V ch through l(TK|l , while the non-negative term 
Q± is linked also to the vector D± through (|17fl . which, 
in turn, contains both the drift of the Gaussian channel 
and the drift created by Bob's additional displacement 



6± 



e(-R) 



Without such a displacement (i.e. 



S± = 0), the teleported state acquires a nonzero drift 
from the channel and the corresponding fidelity -F±(0) 
will depend upon such a drift via a decreasing exponen- 
tial. This does not happen if D ch = 0, or, more generally, 
if g?i — c?3 = g?2 — di — 0. In the general case, the only way 
in which Bob can eliminate the nonzero drift is to choose 
an additional displacement S± which perfectly cancels the 
effects of D ch , i.e., is such that D± = 0. Using (JTSJ, that 
means: 



=Fdi — g?3 and 5 



±di 



(19) 
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With such a displacement the teleportation becomes in- 
dependent from the channel drift and the fidelity 1141) 
becomes 



F± = 



1 



y/det(E ± )' 



(20) 



Suppose now that our input state is a coherent state and 
that our quantum channel has the CM in the standard 
form 



V 



ch 



( a c 0\ 

a d 

c b 

\0 c! b J 



(21) 



then the fidelity (|2C)|> is simply related to the quadrature 
variances: 



F+ 



1 + (Al|) 1 + (AP 



-1/2 



(22) 



where explicitly (AX^_) = a + b ± 2c and (AP|) = 



6 ± 2c'. In particular, if c' 



then (AX\) 



(AP|) > 1 and F± < 1/2 = F c / ass , i.e., quantum tele- 
portation is not possible. Instead if c' — — c then we have 

(AXl) = (AP 2 ) and F± = [l + (AA|)1 . In such a 

case the existence of EPR correlations © is equivalent 
to quantum teleportation, i.e. 



(AX 



< 1 F± > - = F ctess . (23) 



In this sense the Gaussian state $ ch (fli,ftj) < — ► V ch , D ch 
having the CM of Eq. (J2H) with c' = — c is an EPR chan- 
nel generalizing the two-parameter EPR channel studied 
in [lj ■ Note that the left hand side of Eq. (12311 is a 
sufficient condition for bipartite entanglement [K 



IV. TRACING OUT ONE MODE 

We can consider our optomechanical system (Fig. 
as a 3-mode teleportation network where we distill 
a quantum channel between two arbitrary modes i 
and j tracing out the remaining mode k (Fig. 0). 
We denote with — trkipou) the reduced state. 

The corresponding characteristic function is given by 
$W fij) = $(/j,i,(ij,iik = 0), and using Eq. J3J), we 
obtain $ (fe) (/*j; Uj) = cxp[-(/2,, pj)V^(fli, ftj) T ], where: 



( Q l 





H k T k 





\ 





Qi 





-T k 




{-) k T k 





Qj 
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-T k 





Q.7 


/ 



(24) 



that Alice (i) and Bob (j) can arrange the teleportation 
protocol of Section HTD without any additional displace- 
ment 6± and with fidelity F± given by Eq. 1)22(1 for a 
coherent input. It is easy to see from Eq. (|24|) that for 
k = 1 we have F± < 1/2 while for k = 0, 2 we have an 
EPR channel with fidelity = [1 + Q 4 + ± 2T fc ]- 1 . 

In Fig. El we report ^ 0) ,^ 2) versus t' for n = 0, 10 3 
and choosing r = 1 + 2.5 x 10 -7 as in Ref. 0- We 
have F^\2tt) = 1/2 but while f[ 0) < 1/2 Vt', we have 
Fl 0) > 1/2 W (at least up to n = 10 5 , see also Fig.EJ. 
In particular, is approximately constant near 2tt and 

has the expansion Fl 0) - [2 - (*' - 2tt) 2 /2}- 1 . Analyz- 

(2) 

ing Fj_ ; , we see that the quantum channel distilled for 
the Stokes mode (1) and the mirror mode (0) has EPR+ 
correlations before 2w (exploited in [5J) and EPR— corre- 
lations after 2n. The existence of this symmetric region 
is important not only because it allows to perform an 
additional quantum teleportation just after t' = 2n (<^> 
t' = 0) but also because in the time interval delimited by 
the intersections between F^ and F^ (see Fig- these 
latter fidelities are strictly greater than 1/2. This means 
that a region exists after 2ir where a telecloning proto- 
col can be arranged using the Stokes mode (1) as "port" 
and the anti-Stokes (2) and the mirror (0) modes as "re- 
ceivers" [Alice (1) measures (x-,p+) and classically sends 
the result j'_ to Bob (0) and Charlie (2) , who perform the 
displacements ctj — ■> a.j—-f'_,j = 0, 2]. The clones created 
at the receivers are both quantum, but the asymmetry of 

the total channel is such that F^ is only slightly greater 

(2) 

than 1/2 in all the "telecloning interval" , while F_ takes 

its maximum [Fi 2) ] max = (1 + r) 2 /[l + 2r(l + r)] ~ 0.8, 
Vn, and for t' maK = <r/2 + 2tt with <r = cos _1 (2r -2 - 1). 

Since f_^(C- ?/2) = Fl 2) (C ax + ?/2) = (2 + n)~\ 
the telecloning interval becomes narrower by increasing 
the temperature. On the other hand, the behavior of F^ 

(see F^°' tr ^ in Fig.|SJ shows the robustness with respect to 
temperature of the EPR— correlations relative to the dis- 
tilled channel between the two optical modes: the curves 
for n = and h = 10 5 are almost indistinguishable. Such 
EPR— channel, which reduces to a TMS state at t' = ir 
(it is in particular F (Q) (tt) = 1/2 + r/(r 2 + 1) ~ 1), 
enables quantum teleportation even at high tempera- 
tures and points out the present optomechanical system 
as an alternative source of two-mode squeezing. From 

the r-dependence of the maximum values [F_ ] max and 

, we can see that the optimal values, achieved for 
r — > 1, are 4/5 and 1 respectively. 



HETERODYNING ONE MODE 



In other words the distilled channel is a zero-mean Gaus- 
sian state with CM in the standard form l|24|l . It follows 



Consider now our optomechanical system as a 3-mode 
teleportation network where we distill a quantum channel 
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between two arbitrary modes i and j heterodyning the 
remaining mode k and sending the result to Bob through 
a classical channel (Fig.EJ). We denote with a the mea- 
surement result and again with the state of the re- 
duced system involving modes i and j. The correspond- 
ing characteristic function is given by {pi, p,j) — 

N k (a) j dVfc<Q0iG«fc)|a>$(M) with N k (a) = n^ 1 (Qk + 
l/2)exp[(Q fc + 1/2)- 1 \a\ 2 ]. Using Eq. Q, we 
obtain = exp[-(/Zj, Mj)F (fe) Mi) T + 

L> (fc) (^i, /Xj) T ]) where the CM U (fc) is given by (JHJl except 
for the replacements 



U) 



T 2 



Tk. — > Ti. 



TTj 



1/2 



(25) 



(26) 



Ofe + 1/2- 

and the drift is given by 

D (k) = 2i 

Qk + 1/2 

x H^'^T^/'T,) 

[Notice that depends on a and on the i — j order- 
ing too] . It is evident that the distilled channel is still a 
Gaussian state with CM in the standard form but with 
a non-zero drift. It follows that Alice (i) and Bob (j) 
can arrange the teleportation protocol of Section lTTTl with 
a suitable additional displacement 0j_ and with fidelity 
F± given by (|22|l for a coherent input. For k = 1 we 
have again F±^ < 1/2, while for k = 0, 2 we have an EPR 

channel with a modified fidelity F± = [1 + Qi + Qj ± 
2Tk-(T i TT J ) 2 (Q k + l/2)- 1 }- 1 . In such a case (k = 0,2) 

we have S ( ± } = a{T i TT ] )(Q k +l/2)- 1 . In Fig. Hwe report 

versus t' and for n = 0, 1, 10 7 (r = 1 + 2.5 x 1(T 7 ). 



When n — we have Fi 

n(2 



> 1/2 for 7r < t' < 27r and 

F^> > 1/2 for < t' < 7T. However, when n ^ 0, the 
quantum character survives only in small intervals before 
t' = 2ir and after t' = (or equivalently after t' = 2tt). 
These fidelities not only have a temperature-independent 
maximum value [J^ ] max ~ 0.85 greater than that of the 
traced out mode case, but also the "quantum time inter- 
vals" are larger and more robust with respect to tem- 
perature. Moreover, as in the traced out mode case, the 
distilled channel shared by the Stokes mode (1) and the 
mirror mode (0) has both types of EPR correlations and 
therefore one can perform a further quantum teleporta- 
tion just after t' = 2ir (<=>■ t' = 0) in addition to the 
one suggested in 5]. Finally, in Fig. we compare F_ 

versus t' when the vibrational mode is detected (F_ ,het ^) 
with the corresponding fidelity when the same mode is 



distill more EPR correlations between the two optical 
sidebands. In other words, the heterodyne measurement 
allows to improve the efficiency of the present scheme as 
an alternative source of two- mode squeezing • 



VI. CONCLUSIONS 

In conclusion, we have presented an optomechanical 
system as a paradigm of three-mode teleportation net- 
work. We have provided a thorough study of entangle- 
ment and teleportation capabilities of this quantum chan- 
nel. The teleportation fidelities result improved by using 
heterodyne measurement at the remaining mode and us- 
ing the acquired information to distill a finer channel. 

On one hand, our results could be useful to extend 
quantum information processing towards macroscopic 
domain, by using e.g. micro-opto-mechanical-systems 
[l4|. On the other hand, they could be applied as well 
to all o ptic al systems described by the Hamiltonian of 
Eq. O 1I3- 
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FIG. 1: Scheme of the optomechanical system. 



Alice 




Bob 



traced out (F_ 



(0,tr) 



FIG. 2: Three-mode teleportation network. In one case mode 
k is simply traced out while in the other case it is heterodyned 
and the measurement result a sent to Bob. From a zero-mean 
Gaussian state for i,j, k we always distill a bipartite Gaussian 
state. In the first case the quantum channel has zero drift 



for n = 0, 10 . The improvement of and Bob performs the displacement <Xj — > aj — 7±. In the 



the fidelity brought by the additional heterodyne mea- 



surement (F_ 



(0,ftei) 



1 for n = 0, except very close 



to 2m7r) is impressive, meaning that heterodyning the 
mirror mode (and communicating the result) allows to 



second case, the channel has a nonzero drift, known from the 
knowledge of a, and Bob performs a modified displacement 
Qj — > aj — y± — S± with S± depending on the drift (see text). 
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FIG. 3: Traced out mode case. Fidelities F± , F± versus 

t' (around t' = 2-k) for n = and n = 10 3 . We have set 

r = 1 + 2.5 x 10 -7 . The telecloning time interval has been 
marked. 




FIG. 4: Heterodyne detection case. Fidelities versus t 
(around t' = 2k) for n = (a), n = 1 (b) and n = 10 7 (c). 
We have set r = 1 + 2.5 x 10" 7 . 
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